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We have performed multicanonical simulations to study the critical behavior of the two- 
dimensional Ising model with dipole interactions. This study concerns the thermodynamic phase 
transitions in the range of the interaction 5 where the phase characterized by striped configurations 
of width ft = 1 is observed. Controversial results obtained from local update algorithms have 
been reported for this region, including the claimed existence of a second-order phase transition 
line that becomes first order above a tricritical point located somewhere between 5 = 0.85 and 
1. Our analysis relies on the complex partition function zeros obtained with high statistics from 
multicanonical simulations. Finite size scaling relations for the leading partition function zeros yield 
critical exponents v that are clearly consistent with a single second-order phase transition line, thus 
excluding such tricritical point in that region of the phase diagram. This conclusion is further 
supported by analysis of the specific heat and susceptibility of the orientational order parameter. 
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I. INTRODUCTION 

The two-dimensional (2D) Ising model with nearest neighbor ferromagnetic exchange interaction (J > 0) and 
dipolar interaction (g > 0) presents a rich phase diagram because of these competing interactions. This model has 
been the focus of considerable theoretical interest, and the study of its phase diagram has revealed a variety of unusual 
magnetic properties. Moreover, at atomic level, it may give some insight into the interactions that form the striped 
phases observed in a number of ultrathin magnetic films [U [2] as a consequence of the reorientation transition of their 
spins at finite temperatures. The thermodynamic behavior has been investigated by analytical methods and Monte 
Carlo (MC) simulations, aiming at the determination of its critical behavior as a function of the ratio between the 
exchange and the dipolar interaction parameters, 8 = J / 1 g. The Hamiltonian of this model is written as 

<i,j> i<j y 

The variables cr, = ±1 stand for the Ising spins in square lattices L x L and are supposed to be aligned out of plane. 
Here, we have adopted the convention [2j of summing up over all distinct pairs of lattice spins at distances Vij to 
define the dipolar interaction g. The distances are measured in units of lattice. 

Analytical methods include some approximations like spin-wave theory and mean-field, [4l413| but conclusions like 
the fact that the spontaneous magnetization is zero for all temperatures and that the T = configurations present 
patterns classified as regular checkerboards, irregular checkerboards, or stripes of different widths are important. The 
checkerboard pattern corresponds to the formation of alternate magnetic domains represented by black and white 
rectangles. Each of these rectangles contain sites with identical spins and are denoted by < m, n >, where m and n 
stand for lattice units |12j . Regular and irregular checkerboards are defined for m = n and m ^ n, respectively. The 
striped patterns correspond to the formation of magnetic domains displayed in rectangles of size < m, n > but with 
n — > oo. 

Efforts have been made toward a rigorous theoretical proof for the spontaneous formation of these T = configura- 
tions [HI [H] • The formation of such patterns as a consequence of the long-range character of the dipolar interaction has 
been confirmed by MC simulations in different regions of the phase diagram (6, T) . In Fig. [T]we show the phase diagram 
obtained from MC simulations for the 8 range [0,1.9] where the above described ground-state patterns occur. The 
particular case 8 = 0, a pure dipole interaction model, presents a continuous phase transition with critical exponents in 
agreement with the ones in the universality class of the 2D Ising model [HJ[T5]. For < 8 < 0.4152, the model presents 
antiferromagnetic (AF) ground-states characterized by stable regular checkerboard- like spin configurations < 1,1 >. 
Estimates from the specific heat indicate a continuous thermodynamic phase transition associated with the change 
from this AF phase to a phase with broken orientational order, the so-called tetragonal phase [12] , In this phase, the 
magnetic domains lose their common orientation and try to assume the lattice symmetry. The regular checkerboard 
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FIG. 1. Phase diagram: (o) data from Ref. [12 , (♦) this work. Vertical dotted lines represent the phase boundaries through a 
sequence of phases characterized by antiferromagnetic (AF) ground state configurations, the h = 1 and h = 2 striped phases. 
The continuous line corresponds to the expected second-order transition except in the narrow 8 range [0.4152, 0.4403], and the 
dashed line ( ) refers to a first-order one according to Ref. [T2] and [15] . 
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configurations change to irregular checkerboard-like configurations < l,n > in the narrow range 0.4152 < 5 < 0.4403. 
In this 5 range, a thermodynamic first-order phase transition seems to take place |12j . For larger 6 values, the ground 
state changes to spin configurations characterized by magnetic domains displayed in stripes of alternating spins, whose 
stripe width h increases with 6 [HI U3] • Striped configurations of width h = 1 and h = 2 occur for 0.4403 < 5 < 1.2585 
and 1.2585 < S < 2.1724, respectively. Figures [2] and [3] contain these magnetic patterns obtained from our simulations 
for the couplings S — 1.20 and S — 1.30, respectively. In figure [2ja), the low-temperature (T = 0.270) simulation at 
5 = 1.20 presents stripes of width h = 1. Our simulations indicate a transition from the striped to the tetragonal phase 
at T c = 0.311 (figure [2jb)). The tetragonal phase is depicted in Figure [2jc). Figure [3] presents magnetic patterns 
from simulations performed at S = 1.30, a region where stripes of width h = 2 occur. 

In addition to the striped and tetragonal phases, a new domain in the phase diagram has been reported for 5 = 2, 
the so-called nematic phase in analogy with liquid crystals. In the nematic phase the system still keeps its long-range 
orientational order but loses the spatial order exponentially. This new domain has been studied by MC simulations, 
and it has been found between the striped and the tetragonal phases. In this case, we would have two thermodynamic 
phase transitions: stripe-nematic and nematic-tetragonal transitions. This new phase is located in a region of the 
(S,T) plane that gives origin to a bifurcation of the line that separates the h = 2 and h = 3 phases p~6| I17j. 

A convincing determination of the thermodynamic phase-transition order is still lacking even for such small h 
values. In fact, controversial results about the order of the thermodynamic phase transition as a function of 5 have 
been reported in the literature. In particular, some MC results concerning square lattices for 5 between 0.2 and 2 
exhibit a phase diagram with a second-order transition line ^] for the thermodynamic transition between the 
ordered phases and the tetragonal one. On the other hand, the transition line appears to be first order for a 6 range 
corresponding to h > 1 [T2], with the remark that for 6 = 0.85 a second-order phase transition takes place with 
exponents dv — 2.0 ± 0.1, a = 0.09 ± 0.07, and 7 = 1.75 ± 0.05 at the critical temperature T c = 0.41 [15 . As usual, 
v, a, and 7 are the correlation length, specific heat, and susceptibility exponents, respectively. For the interaction 
8 = 1, it seems to be first order at T c = 0.40 [12] and to present only a weak first-order character at T c = 0.404 
|16| . The above cited results would lead to the existence of two critical lines separated by a tricritical point for 6 
somewhere between 0.85 and 1. [7] 

The controversial results are a consequence of the dipolar term, which produces large autocorrelations in MC time 
series obtained with local update algorithms |17H19j . Moreover, simulations have also been hampered because large 
lattice size simulations are very CPU time consuming due to this term, frustrating any convincing finite size scaling 




FIG. 2. Spin configurations for 5 = 1.20 and L = 72. (a) striped phase: T = 0.270, E/N = -0.4638, O hv = 0.9869; (b) 
transition temperature: T c = 0.311, E/N = -0.4096, O hv = 0.5039; (c) tetragonal phase: T = 0.350, E/N = -0.3539, 
O hv = 0.0476. 




FIG. 3. Spin configurations for 8 = 1.30 and L = 72. (a) striped phase: T = 0.290, E/N = -0.5058, O hv = 0.9831; (b) 
transition temperature: T c = 0.329, E/N = -0.4397, O hv = 0.5173; (c) tetragonal phase: T = 0.380, E/N = -0.3821, 
O hv = 0.0370. 



4 



(FSS) analysis. In this paper, we perform extensive multicanonical simulations for determination of the character of the 
thermodynamic transition from the h = 1 phase, and to address the existence of a tricritical point. The multicanonical 
algorithm (MUCA) generates a ID random walk in the energy space, diminishing the problem of overcoming free 
energy barriers. We carry out a comprehensive analysis of the character of the phase transition for values of S from 
0.85 up to 1.30 by means of the complex partition function zeros [2H1 HI]- Partition function zeros analysis in the 
complex temperature plane has been successfully applied to spin models [22-24 , lattice gauge theories [25J ES] , and 
protein models [571125]. This procedure has allowed us to explore critical aspects by means of FSS relations for the 
first complex zero, leading to a precise characterization of the phase transition line. The conclusions based on the 
partition function zeros are further supported by analysis of the specific heat and susceptibility of the orientational 
order parameter. Analysis of these thermodynamic quantities allows us to calculate the critical exponents ajv and 
7/V. It is well known that the renormalization-group fixed point picture for d-dimcnsional systems in the L d block 
geometry characterizes a first-order phase transition by the particular value of the critical exponent dv — 1 |291 130j . 
This, in turn, gives a — 1 and 7 = 1 for a first-order phase transition, which produce the expected dependence of the 
thermodynamic quantities on the volume L d and has been supported in a number of Monte Carlo studies [24l |3T] . 

Our results rely on data collected from lattice sizes up to I = 72. We report precise estimates for the infinite 
volume critical temperatures and critical exponents v, a/v, and j/u for values of 6 from 0.85 up to 1.20. We have 
included the interaction <5 = 1.30 in the h = 2 phase for comparative purposes. In Sec. [TTJ we briefly review the main 
aspects of MUCA, and the protocol devised for updating of the multicanonical parameters. In Sec. |III[ results from 
FSS relations for the first complex zero, specific heat, and susceptibility are compiled, to produce the estimates for 
the critical exponents. The final Sec. IV presents a summary and our main conclusions. 



II. MULTICANONICAL SIMULATIONS 



The multicanonical algorithm |321 133] , like other generalized algorithms |34j , significantly improves the sampling 
of configurations. This algorithm assigns a weight w mu (E) ~ l/n(E), where n(E) is the density of states and 
E = H({<ii}) is the energy of a state given by the spin configuration {o^}, with i — 1, ...,L 2 , as defined in Eq. |l]). 
Therefore, the multicanonical method is expected to produce flat energy histograms H mu (E) cx n(E)w rnu (E) under 
the following probability condition 



p(E -» E') = 



W m u(E') 

' w mu (E) 



(2) 



for sufficiently long simulation times. 

The multicanonical weight w mu (E) is a priori unknown. A numerical estimate of w mu (E) is usually obtained by 
considering the Boltzmann entropy S(E) = \nn{E) (ks — 1), and the following parameterization for the entropy 
S(E) = b(E)E — a(E), where a(E) and b(E) are called multicanonical parameters. Hence, the multicanonical weight 
is given by w mu (E) — exp[—b(E)E + a(E)], with the parameter a(E) related to a multicanonical free energy and b(E) 
related to the inverse of the microcanonical temperature. 

The implementation of MUCA requires energy discretization. An integer label m is introduced to facilitate our 
histogramming of energy data. This label defines energy bins of size e, E m = Eq + me, with m = 0, • • • , M. All 
the energies in the interval [E m , E m+ i[ are in the mth energy bin and contribute to the histogram H mu (E m ). The 
constant Eq is defined as a reference energy just below the ground-state energy. We have verified that e = 1 is a 
convenient discretization. 

The parameters a(E) and b{E) are estimated from N r recursion steps. Each step updates the multicanonical 
parameters through the following equations [32] . 

a n {E m _ x ) = a n (E m ) + [b n (E m ^) - b n (E m )]E m , 

b n (E m ) = b n -\E m ) + [\nH^(E m+1 ) - \nH^{E m )]/s , (3) 

where n (n — 1, • • ■ , N r ) labels the recursion steps and N r amounts to how long this update procedure is enforced 
in order to obtain reliable estimates for w mu (E). The choice 7J^ u (i? m ) = m&x[ho, i?^ lt( (£ , m )], with < ho < 1 
for all discretized energies, is a technical choice to avoid H(E m ) = [32] . It is convenient to compute the above 
recurrence relations with the initial conditions a°(E m ) = and small values for b°(E m ) if the simulation uses a hot- 
start initialization. The nth recursion step needs the calculation of H mu (E m ) from the previous weight {a™~ 1 , 
obtained with n s MC sweeps. Usually, the number N r is defined a posteriori when the multicanonical parameters 
present some convergence. 

To determine the multicanonical parameters, we have devised the following protocol. Each recursion step is 
implemented after collection of H mu data by sampling configurations between two extremal energies E*_ and E*^, 
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with E*_ < E+. A round trip is defined as the number of sweeps necessary to go from configurations with the lowest 
reference energy E*_ to the ones with a fixed high energy E* + and back. A round-trip walk may also start at any 
energy between E*_ and £7*. The multicanonical update procedure Eq. (|3j) is performed with a variable number of 
MC sweeps necessary for the attainment of three of such round trips. This number of round trips is chosen to ensure 
samplings across the energy landscape in a reasonable simulation time. To avoid too long simulation time to achieve 
the next (n + l)th multicanonical recursion, a fixed number of sweeps n s (n) is set as the limiting number of MC 
updates. Thus, new multicanonical parameters are obtained as soon as one of the following conditions is observed: a) 
three round trips or b) a number of MC sweeps greater than three times the average number of MC sweeps counted 
in the previous multicanonical simulations, 



n-l 



•(») = — ,E^w- ( 4 ) 

71—1 



1=1 

After each multicanonical update, E*_ is replaced with the minimum energy among the sampled energies in the 
previous simulation. This establishes a new (and larger) energy interval where new round trips must occur. This 
protocol helps us to keep a reasonable number of tunneling events even for large lattice sizes at the price of longer 
CPU times. A further improvement of the multicanonical weight w^ u is achieved with an extra MUCA update, which 
consists of timc MC sweeps necessary for the performance of 20 round trips. Table I lists only the number of sweeps 
fiMC a s a function of the lattice size L for different interactions S that are necessary for the accomplishment of this 
final update. With this final estimate of the multicanonical weight w mu {E), we proceed to data production. Our 
data production amounts to 16 independent energy time series, each one produced with tlmc sweeps. Thus, the data 
analysis for the smallest lattice size L — 12 and 5 — 0.89 relied on ~ 1.8 x 10 5 measurements, while in the case of the 
largest lattice size L = 72 and 6 = 1.30 it amounted to ~ 1.62 x 10 s measurements. We can anticipate that the large 
number of measurements for 6 = 1.30, compared with the smaller 5 values, is related to the effort of overcoming the 
free energy barrier as a consequence of a first-order phase transition at this interaction. 

We have carried out simulations with periodic boundary conditions to minimize border effects. Thus, all distances 
Tij must include sites in the infinitely replicated simulation box in both directions. This boundary condition adds an 
infinite sum over all images of the simulation box because of the dipole term in the Hamiltonian. The infinite sum 
was computed by means of the Ewald summation technique. This technique splits the infinite sum over all images of 
the system into two quickly converging sums, namely the direct sum, which is evaluated in the real space, and the 
reciprocal sum, carried out in the reciprocal space, as well as a self- interaction correction term [33 [36] . We set the 
Ewald parameter a to 3.5 in all the simulations. This parameter determines the rate of convergence between the two 
sums. 

An important consequence of MUCA data production is the estimation of canonical averages of thermodynamic 
quantities A over a wide range of temperatures T — 1//3 by using the reweighting technique [37j 



E fc A k [w^ u {E k )]-^{-PE k ) 
Efe[^m«(^fe)]" 1 exp(-/3^ fe 



A{fi) = ^ k \j mu " /J ^ v - Kl . (5) 



This contrasts with the Metropolis algorithm, where the reweighting is restricted to a very narrow temperature range 
around the fixed MC simulation temperature. After reliable estimates for the MUCA weight, one can evaluate the 
density of states 

n(E) = H mu (E)w^ u (E), (6) 
from which one can construct the partition function 

Z((3) = J2n(E)u E , (7) 

E 

where u — eT^ . The complex solutions in u, {Re(u), Im(u)}, describe the critical behavior of the system. These 
solutions correspond to the so-called Fisher zeros [20] HI] • 

III. RESULTS 

A. Partition function zeros 

Let us consider the complex zeros of Eq. Q ordered according to their increasing imaginary part. For a sufficiently 
large lattice size L, the leading partition function zero u\(L) can be used to obtain the critical exponent v through 
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the FSS relation [3T], 

u?(L) =u c + AL- 1 ' v [l + 0{L y )\ , y<0. (8) 

This relation shows that the distance from the closest zero to the infinite lattice critical point u c — e~@ c on Re(it) 
scales with the lattice size. If we disconsider finite size corrections, the exponent v can be obtained from the linear 
regression 

- \k\u\(L)-u c \ = -hx(L)+a . (9) 
v 

Since the exact critical temperature is unknown, and because the real part of u presents weaker dependence on L 
as compared to the imaginary part of u, it is usual to replace |it° — u c \ with its imaginary part, so as to avoid a 
multiparameter fit. 

With the discretization e, Eq. (|7|) becomes a polynomial in u and it can be solved with MATHEMATICA for 
L < 32. Larger lattices present huge numbers for the density of states, which makes the scan method in the complex 
temperature plane the only way of obtaining complex zeros |38j . The leading complex zeros are presented in Tables 
II, III, and IV as a function of L for different S values. 

Now, considering the real part of those zeros, Re^^i)] = — l/21n{[Reu°(L)] 2 + [Imii5(L)] 2 }, one can estimate 
the critical temperatures through the following FSS fit [31] : 

Re \ft{L)] = /3 L ° + bL- 1 ^ . (10) 

This fit yields the critical temperatures T~ displayed in Table V, where we have included the exponents dv. 

The quality of the linear fits can be stated in terms of the goodness-of-fit of the model [39]. The goodness-of-fit Q, 
(0 < Q < 1) is related to % 2 and, as a general rule, if Q is larger than 0.1, then the fit is believable. The linear fits 
for evaluation of v present Q as large as 0.98 for 5 < 1.10, Q = 0.70 for 5 = 1.20, and a very small Q value ~ 10 -10 
for S = 1.30. Figure [4] illustrates our linear fits for S > 0.97. The data fit nicely, confirming the linear dependence on 
ln(L). However, the very small Q value for S = 1.30 seems to be a consequence of high statistical precision for these 
zeros (see Fig. [4]), which reveals the presence of some systematic bias. It is known that corrections to FSS relations 
give a better fit for first-order phase transitions [31 j - However, this would require larger lattice sizes for the attainment 
of reliable estimates from a multiparameter fit. The asymptotic behavior of Ti(L) is determined with high Q values 
for all 6 interactions. Data collected in the third column of Table V shows a consistent trend toward dv = 1 as we 
move on the critical line in the direction of higher 5. The value dv = 1 is only reached for the interaction in the h = 2 
phase. Thus, these results clearly exclude first-order phase transitions from the h — 1 phase. 

Results for are depicted in Fig. [ljwith the symbol (♦). In this figure we also show the values obtained from 
Ref. [15] and, in particular, we note that the values for 5 = 0.85 and 1.3 are surprisingly good as compared to Tj?, 
since they are obtained from a single lattice size L = 48. 
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FIG. 4. Finite-size scaling fits of the leading complex zeros for some S couplings. 
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B. Specific heat and susceptibility 



To further characterize the order of the phase transitions, we have studied the specific heat, 

1 



C V {T) 



T 2 N 



{{E*)-{EY), 



and the susceptibility 

X(0 hv ) = N «OL) - (Ohv) 2 ) 
associated with the orientational order parameter [40] . 



Ohv 



(11) 



(12) 



(13) 



over a (continuous) range of temperatures by reweighting MC data according to Eq. ([5]). The quantities and n v 
are the number of horizontal and vertical bonds of the nearest neighbor antiparallel spins, respectively. This order 
parameter is +1 in the striped ground state, and it vanishes at high temperatures where orientational symmetry of 
the striped domain is broken. A very common way of obtaining the critical exponents is through the FSS relations 
for the maximum of the specific heat 



and for the maximum of the susceptibility, 



a|max(T c (L),L)aL a /" 



Xmax(T c (L),L) (XL?/* 



(14) 



(15) 



where T C {L) is the finite size critical point. Again, an FSS relation like Eq. (10 1 is applied to the temperatures 
T C (L), where the maxima of C V (T,L) and x(T, L) occur, to yield the infinite volume critical temperature T^ v and 
T c x , respectively. Table V summarizes these temperatures and the critical exponents for C v and \- The temperatures 
If" and T* are then evaluated with v obtained from the hyperscaling relation a = 2 — dv, with data displayed in 
the 5th column of table V. The goodness-of-fit of the linear fit for C v is about 0.5 for S < 1.20. Again, it decreases to 
a very small value Q ~ 10 -5 for 5 = 1.30. The linear fit of \ presents Q ~ 0.8 for 5 < 1.20 and also decreases to 10~ 5 
for S = 1.30. 

The critical exponents ajv in the 5th column (Table V) clearly exclude any possibility of a first-order phase 
transition from the h = 1 phase, while it strongly indicates this possibility at 5 — 1.30. The statistical error bar 
exclude the value ajv — 2 at <5 = 1.30, but the small Q value may indicate the presence of systematic bias. The 
results from the susceptibility are less prompt to make satisfactory claims about the order of the phase transition 
only at 5 = 1.2 and 1.3. Again, those results may be due to the missed corrections to the FSS relation, as expected 
at first-order phase transitions. Figures [5] and [6] display the FSS plots for C V (T C (L),L)) and x(T c {L),L)) for 6 = 1.20 
and 1.30, respectively. These figures have helped us observe how satisfactory the FSS are. 



IV. SUMMARY AND CONCLUSIONS 



We have performed analysis of the complex partition function zeros from multicanonical simulations. The sampling 
with this algorithm is known to be efficient when it comes overcoming the free-energy barrier problem in simulations 
of complex systems as compared to the usual local update algorithms. MUCA is based on a non-Boltzmann weight 
factor and performs a free one-dimensional random walk in the energy space. A protocol has been devised for 
the determination of the multicanonical weight factor by ensuring that enough measurements in the energy space 
are obtained. Moreover, by keeping such control over the number of round trips between the low and high-energy 
configurations, we were able to determine the number of sweeps that is necessary for exploration of the energy space 
even for large lattice sizes. 

By using FSS relations involving the partition-function zeros obtained with high statistics, precise estimates for 
the infinite- volume critical temperatures and critical exponents v were found for interactions S corresponding to the 
h = 1 phase. We also included an interaction (5 = 1.30) that produces stripes of width h = 2 for comparative 
purposes. Analysis of the specific heat and susceptibility of the order-disordered parameter Ohv gives further support 
for a second-order transition critical line from the h = 1 phase. Infinite volume critical temperatures obtained from 
ul, maximum of C v and Xi are in full agreement, which help us draw a reliable part of the phase diagram (5, T). In 
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FIG. 5. Finite-size scaling plots for the (a) specific heat and (b) susceptibility as a function of the temperature for 5 = 1.20. 
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FIG. 6. Finite-size scaling plots for the (a) specific heat and (b) susceptibility as a function of the temperature for 5 = 1.30. 



9 



conclusion, the study conducted with many 5 interactions strongly indicates the existence of a second-order critical 
line between the high-temperature tetragonal phase and the low-temperature ordered phase characterized by h = 1. 
The first-order character is found in our study only for the interaction 5 = 1.30. This suggests that the second-order 
critical line ends at 6 = 1.2585, and that it becomes first-order beyond this point. We assume that both transition 
lines are separated by a tricritical point at S = 1.2585, because this point produces a line separating the h = 1 and 
h = 2 phases, and it bifurcates for generation of the tetragonal phase. The precise temperature where this bifurcation 
happens has not been evaluated in the literature to the best of our knowledge. 
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5 0.85 


0.89 


0.91 


0.93 


0.95 


0.97 


1.00 


1.10 1.20 


1.30 


12 


15 018 


11 632 


12 032 


14 026 


10186 


16 893 


12 279 


24216 51171 


73 362 


20 


65 278 


47 290 


52 098 


46 969 


50127 


78 380 


63 333 


91158 226 934 


269 598 


32 


181069 


188 620 


149 412 


105 008 


158 214 


213 036 


179 126 


236 436 959 554 


1354 207 


48 


401 195 


483 287 


380 451 


645 854 


503 469 


484 012 


582 030 


879 740 2 285 886 


3 256 588 


64 


1 028 097 


824 022 


938 475 


751 156 1 104 830 1 300 445 1 296 208 2 281 061 5 579 023 


8159 701 


72 


1 283 737 1 345 861 1 421 650 1 332 818 1 709 722 1 590 101 1 697 157 1 975 543 5 682 748 10 171 715 



Table I. Number of MC sweeps umc as a function of the lattice size L for different interactions 8. These numbers correspond 
to 20 round trips observed in the final MUCA update. 





8 = 0.89 




8 = 0.91 




8 = 0.93 




L 


Re(u?) 


Im(u?) 


Re(w?) 


Im(u?) 


Re(u?) 


Im(u?) 


12 


0.0913(37) 


0.0401(28) 


0.0923(36) 


0.0381(29) 


0.0907(34) 


0.0377(21) 


20 


0.0911(22) 


0.0222(14) 


0.0909(24) 


0.02117(95) 


0.0899(14) 


0.0205(13) 


32 


0.0908(12) 


0.01330(55) 


0.08967(96) 


0.01300(90) 


0.0886(11) 


0.01237(62) 


48 


0.09035(67) 


0.00842(57) 


0.08869(99) 


0.00807(51) 


0.08776(49) 


0.00770(50) 


64 


0.08932(99) 


0.00615(39) 


0.08864(64) 


0.00598(32) 


0.08729(54) 


0.00568(39) 


72 


0.08970(60) 


0.00556(39) 


0.08857(85) 


0.00522(41) 


0.08711(64) 


0.00502(31) 



Table II. Complex partition function zeros for 8 = 0.89, 0.91, and 0.93. 





8 = 0.95 




8 = 0.97 




8 = 1.00 




L 


Re(it?) 


Im(w?) 


Re(it?) 


Im(it?) 


Re(u?) 


Im(it?) 


12 


0.0895(31) 


0.0358(25) 


0.0871(21) 


0.0339(16) 


0.0844(18) 


0.0312(16) 


20 


0.0886(14) 


0.0194(10) 


0.0858(11) 


0.01841(71) 


0.08298(86) 


0.01637(41) 


32 


0.08708(82) 


0.01157(76) 


0.08497(55) 


0.01053(47) 


0.08162(62) 


0.00916(42) 


48 


0.08634(51) 


0.00728(34) 


0.08446(42) 


0.00677(41) 


0.08124(43) 


0.00577(23) 


64 


0.08580(27) 


0.00513(17) 


0.08408(28) 


0.00479(28) 


0.08096(17) 


0.00406(25) 


72 


0.08583(31) 


0.00460(25) 


0.08391(26) 


0.00408(21) 


0.08068(15) 


0.00353(14) 



Table III. Complex partition function zeros for 8 = 0.95,0.97, and 1.00. 
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5 = 1.20 5=1.30 
Re(w?) Im(uS) Re(u?) Im(w?) 



6 = 1.10 

L Re(w?) Im( M ?) 

12 0.0666(14) 0.02223(87) 

20 0.06579(63) 0.01049(34) 

32 0.06524(48) 0.00544(21) 

48 0.06504(21) 0.00312(12) 

64 0.06488(10) 0.002070(83) 

72 0.06476(11) 0.0017227(77) 



0.0365(14) 0.01504(51) 
0.03968(82) 0.00687(35) 
0.04006(22) 0.00326(17) 
0.040060(91) 0.001584(91) 
0.039935(59) 0.000974(49) 
0.039934(39) 0.000800(47) 



0.0515(10) 0.01329(29) 
0.05091(45) 0.004203(34) 
0.04876(14) 0.001745(15) 
0.048100(66) 0.0007782(50) 
0.047899(44) 0.0004346(18) 
0.047868(29) 0.0003433(11) 



Table IV. Complex partition function zeros for S — 1.10, 1.20, and 1.30. 



s 


rjiQ 

C 


dv 


■*■ c 


a/v 


T* 


-y/v 


0.85 


0.41189(53) 


1.837(76) 


0.41240(48) 


0.344(16) 


0.41200(51) 


1.519(19) 


0.89 


0.41168(53) 


1.807(70) 


0.41104(62) 


0.364(20) 


0.41100(48) 


1.531(27) 


0.91 


0.40887(50) 


1.817(68) 


0.40992(16) 


0.375(19) 


0.40964(17) 


1.538(22) 


0.93 


0.40681(19) 


1.779(61) 


0.40685(46) 


0.399(20) 


0.40682(45) 


1.561(24) 


0.95 


0.40435(17) 


1.741(53) 


0.40475(12) 


0.424(20) 


0.40475(18) 


1.552(24) 


0.97 


0.40108(40) 


1.706(46) 


0.40124(45) 


0.461(19) 


0.40130(31) 


1.575(20) 


1.00 


0.39499(37) 


1.659(37) 


0.39521(33) 


0.522(17) 


0.39527(29) 


1.590(23) 


1.10 


0.36429(29) 


1.415(25) 


0.36441(19) 


0.888(21) 


0.36456(15) 


1.736(21) 


1.20 


0.31102(32) 


1.223(21) 


0.31126(65) 


1.496(28) 


0.31073(40) 


1.987(29) 


1.30 


0.32929(72) 


1.0093(28) 


0.32892(15) 


2.0183(66) 


0.32885(14) 


2.3193(82) 



TABLE V. Critical temperatures and critical exponents from complex partition function zeros, specific heat C v , and 
susceptibility x{Ohv)- 



